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A proposal is made for reducing the solution of the iV-particle Lippmann-Schwinger equation to 
that of smaller sets of particles. This consists of first writing the iV-particle equation in terms of all 
possible iV/2-particle Lippmann-Schwinger equations. (If N is odd this needs a minor modification.) 
The second step requires a decoupling of the resolvents for the fewer particle systems so that each can 
be solved separately. This generalization of the Faddeev approach deals only with connected kernels 
, and the homogeneous solution reproduces the iV-particle Schrodinger equation. For four particles 

■ the proposed method involves only a 3 x 3 matrix whereas other approaches typically require the 

. solution of at least a 7 x 7 matrix equation. 
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I. INTRODUCTION 



^ ' There has been an increased interest in many-body scattering dynamics due to the rapid growth of computer power 
|fflj. The fundamental method of formal scattering theory is the Lippmann-Schwinger equation pi. This equation is 
equivalent to the Schrodinger equation for a two-body problem with an appropriate scattering boundary condition. 
, However, it is well recognized that there are problems when solving the Lippmann-Schwinger equation for more 
than two particles. In particular for three particles, but the same problem occurs when there are more particles, a 
perturbation expansion in powers of the potential, or better in terms of pair particle T-matrices, allows a third particle 
■ to be a spectator and a dangerous delta function to arise in the momentum representation of the Lippmann-Schwinger 
equation || . Faddeev j| rearranged the equation so as to eliminate the dangerous delta function and in the process he 
was able to show that the kernel of the rewritten integral equation was compact, essentially by showing that it became 
connected after one iteration. Lovelace H and later, Alt, Grassberger and Sandhas (AGS) [|| gave two different ways 
of rewriting the Faddeev equations in terms of scattering amplitudes for various collision processes. The purpose of 



this work is to propose a way of generalizing the Faddeev equation to an arbitrary number of particles. 

There has been a large effort in extending Faddeev's work to an A^-particle system, in particular to the four-particle 
system, both for theoretical interest and for practical use in nuclear and chemical reactions [j3p|- |2l]]2^ , |24|| . Most 
work puts the emphasis on deriving eq uations with a connectable kernel, namely, a kernel that is connected after a 
finite number of iterations. Glockle JT3] showed that a set of generalized Lippmann-Schwinger equations has a unique 
solution for the three-particle system. He also formulated similar equations for a four-particle scattering system 
p5j. Special formulations were given by Weinberg ||| and Rosenberg ffl, the former emphasized the quasiparticle 
picture of an interacting A^-particle system, the latter provided a set of equations that do not explicitly depend on the 
scattering potentials, a feature that may be useful for rclativistic scattering J26|. A channel array approach was taken 
by Baer et al. , to achieve a connected kernel for a three-body system which was later generalized to an iV-particle 
scattering system by Tobocman [^0|. Yakubovsky |l| presented a formal generalization of the Faddeev equations to 
A^-particle scattering by appropriately organizing pair-particle interactions. Yakubovsky 's equations jljj, like those 
of Faddeev, have considerable advantage in dealing with the bound state problem since they are equivalent to the 
Schrodinger equation, thus there being less chance of spurious solutions p^ , p8| . Unfortunately, the complexity and the 
lack of a simple method of iterating the Yakubovsky equations hinder their practical usage in physical and chemical 
problems. In a later treatment, Narodetsky and Yakubovsky |]l4]j proposed an alternative two-cluster approach and 
obtained a new set of equations for the A^-particle scattering. The equivalence of these two-cluster equations to the 
Schrodinger equation has not be proved. Sloan made an important contribution to solving the four-particle problem 
by allowing only two-body channels in the coupled equations |l5|] for scattering amplitudes. Bencze obtained 
generalized Lovelace equations and generalized AGS equations by resuming the two-cluster equations obtained by 
Narodetsky and Yakubovsky and showed that these equations are equivalent to those of Sloan when N = 4. Redish 
p7[ generalized Sloan's equations to an arbitrary number of particles and demonstrated that his results are identical 
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to those of Bencze. The Sloan-Bencze-Redish equations are connected after one iteration provided that the kernels of 
the subpartitions are connected. 

These previously presented methods appear to emphasize the pair particle transition operator and expand everything 
in terms of it. This approach of building up everything from pair particle properties involves combining the pair 
transition operators in more and more elaborate ways, the details of which may often be very tedious. In fact, there 
are so many terms that they are very cumbersome to implement. Even the comparison of the numerical calculations 
using different methods appears to raise questions. What is proposed here might be thought of as an approach of 
pulling the problem apart into smaller pieces. Thus the V-particle system (if N is even) is divided into all possible 
combinations of iV/2-particle subsystems with the coupling of the subsystems carried out in a manner parallel to the 
method of Faddeev. This is repeated for each TV/ 2-particle subsystem. Eventually, after successive divisions, each 
subsystem consists of two particles, whose collisional effects are described by a pair transition operator. (A minor 
modification is required if N is odd.) This gives a structure as to how the iV-particle transition operator is dependent 
on the transition operator of fewer and fewer particles. The calculation of the V-particle transition operator then 
involves the reverse procedure of successively putting together the transition operators of fewer particles to eventually 
get an expression for the full transition operator. In going up the chain of partitionings, each step leads to an integral 
equation of Faddeev type which is connected after one iteration. In this way a set of equations is obtained in which all 
particles are connected, but the connectedness is organized in sets within sets. A necessary technicality for carrying 
out this approach is the need to renormalize the potential at each partitioning. Thus the transition operators for 
smaller sets of particles involve a potential which is generally smaller in magnitude than the true potential. It is 
believed that the approach presented here has a different basis of approach than those previously presented, has fewer 
terms and hopefully leads to a more efficient numerical procedure. 

This paper is divided into five sections. Section II is devoted to the partitioning of A-particle systems, with emphasis 
on partitions having two clusters with an equal number of particles. The notation is similar to that of Redish M] 



and Yakubovsky 
III to generalize t 



|13|. The distinction between a channel and a partition is stressed. A scheme is developed in Section 
;he Faddeev equations to a set of equations for A-particles. The derivation given for the generalized 
Faddeev equations closely parallels Faddeev's original derivation. This emphasizes the simplicity of the new equations. 
The generalized Faddeev equations for the resolvent operator and the decomposition of the scattering wave functions 
are given in Section IV. Section V discusses how the proposed procedure would be applied to both 4- and 8-particle 
scattering. A general discussion in Section VI ends the paper. 



II. PROPERTIES OF EQUAL PARTITIONING 

Consider a system of N labeled particles (1, 2, 3, • • • , N). There are many topologically distinct ways of partitioning 
the labeled N particles into sets of clusters of particles. Each particular way of dividing up the N particles will be 
called a partition, labeled by A, B, C, etc. The detailed presentation given here partitions the N particles into two 
clusters of nearly equal size (if N is odd, one cluster is to have one more particle than the other), but alternate 
partitioning schemes are also possible, as mentioned in the Discussion. All partitionings of the same size are to be 
considered, which differ by particle composition in each cluster, to give a set of partitions C. The objective of this 
section is to define the potentials and hamiltonians appropriate for such partitions. As used in this work, a partition 
only becomes a (rearrangement) channel when the clusters composing the partition are bound, thus distinguishing 
between the mathematical method of solving the iV-particle problem and the physical notion of what asymptotic 
states arise in a scattering process. 

The V-particle hamiltonian H = Ho + V consists of a kinetic energy operator Ho and an assumed pairwise additive 
potential operator 

V = £Vy, (1) 

where the sum is over all possible ordered pairs, with i,j € {1, 2, 3, • • • , N}. The total hamiltonian is assumed to 
be self-adjoint and bounded below, and to have spectra that might include a discrete set, associated with the bound 
eigenstates, as well as the continuum for the scattering states. The objective is to determine the solutions for the 
Schrodinger equation associated with the total hamiltonian H, 

H | *) = £ | (2) 

In the present treatment, the N particles are partitioned into two clusters of equal size if at all possible. That is, 
if N is even, there are N/2 particles in each, whereas if N is odd, one cluster in a partition has (N + l)/2 particles 
and the other (N — l)/2. It follows that the number of partitions in this set C of partitions is given by 
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for N even; 

I ((JV-D/2) for N odd - 

The partition potentials are to be chosen so that 

E V ^ V - (4) 

Thus the potential Vp of partition C is defined as 

V C = %^ V jk , (5) 

j<k 

i.hec 

involving the sum over those pairs of particles jk that appear in the same cluster in the partition C . Here vn is a 
renormalization factor and obviously, W3 = 1. For N > 3, vn is chosen as 



)/[(T)Gvp 



= t^pgbS! iV even; 

^ (?)/ {((, -V) [( (J % 1)/2 ) + ( (JV+ 2 1)/2 )] } - « N+ »'™»- 3) ' 2)l * odd, 

so that the total number of potential terms on both sides of Eq. (||) are equal. The special case of V4, = 1 implies 
that no renormalization is needed for N = 4, whereas D5 = 1/4 and «6 — 1/4 demonstrate how the potential for a 
particular pair of particles gets distributed among the many partitions in the set C when N is large. Note that other 
types of two cluster partitions can also be selected (or included), which however, leads to a larger number of coupled 
equations for a given iV-particle system. 

A partition hamiltonian is defined as the kinetic energy operator Ho plus the mutual potential interactions between 
those particles in each of its clusters, namely the partition potential, thus 

H C = H +V C . (7) 
The residual interaction, V c , describing how the clusters in partition C interact with each other, is defined as 

V C =V-V C . (8) 
Thus, for each partition C the total hamiltonian H can be written as 

H = H C + V C . (9) 
The resolvent operator R(z) for the total hamiltonian is defined as 

R~, do) 

where z is a complex parameter, for scattering theory having a small positive imaginary part. In most equations in 
this work the dependence on z will be implicitly assumed. The free resolvent operator is defined as 

Ro = —^-rr (11) 

z - H 

while the resolvent operator for partition C is 

Rc = (12) 
z — tie 

The object of this paper is to express the total transition operator T, as determined by the Lippmann-Schwinger 
equation 

T = V + VR T, (13) 
and the Schrodinger equation (||) in terms of the clustering discussed above. 
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III. GENERALIZED FADDEEV EQUATIONS FOR T 



In this section a two-step scheme is proposed for the decomposition of the total transition operator T in such a 
manner that its kernel is connected after one iteration. The first step is to decompose T according to the selected set 
of two-cluster partitions C, see Sec. 0. Conceptually, this is exactly what Faddeev Q did for the three-body system. 
However for a system of more than three particles, the resulting expression for a partition transition operator Tc 
is governed by a combination of the transition operators for the two clusters. The second step is to identify how 
this combination is to be carried out, compare Ref. fl5|| . Finally it is recognized that the separate cluster transition 
operators are exactly similar to the original iV-particle transition operator but with significantly fewer particles, so 
the whole procedure can be repeated for each cluster. Care must be exercised when reconstructing the partition 
transition operators from the separate cluster transition operators so that the energies appearing in the resolvents are 
consistent with the total energy of the system. 

A. The iV-Particle Transition Operator 

The full ./V-particle transition operator has the equivalent forms 

T = V + VRV 
= V + TR V 

= V + VR T. (14) 

Moreover, it can be decomposed according to 

T = E v c + E v c R o T 



cec cec 



(15) 



cec 

where 

T C = V C +V C R T. (16) 

Subtracting VcRoT^ from both sides of Eq. (|l~6| ) gives 

(1-V C R )T C =V C +V C R (T-T C '), (17) 

which, on multiplying by (1 — VcRo)" 1 gives 

T c -(l-V c R )" 1 Vc 

+ (1 - VcRq)- 1 V C R (T - T c ) 

= T c + ^l c , B T c RoT B (18) 

Bee 

Here 8c, b = 1 — <5c,s and the (2-cluster) partition transition operator Tc is 

T C = (1-V C R )- 1 V C 

= V c + V C R T C = V c + T C R V C 

-V C +V C R C V C . (19) 

Therefore the whole iV-particle transition operator T can be split into A/jv Faddeev type transition operators T c 
associated with the selected set C of two-cluster partitions. It is interesting to note that Eq. (|lj) has the precise 
structure of the Faddeev equation but is valid for an arbitrary number of particles and with A/jv components. Just 
like the original Faddeev equation, the kernel in Eq. (fL8|) is connected after one iteration provided the scattering 
kernels of the subsystems are already connected. 

As in the other treatments fl5T~|l7j| of the ./V-particle system, 2-cluster transition operators Tc may be expressed in 
terms of unconnected subsystems. Thus it is necessary to ensure that each partition transition operator is connected 
in order for the total transition operator to have a connected kernel. This can be accomplished by assuring that each 
cluster transition operator is connected. But first, it is necessary to know how to express the partition transition 
operators in terms of the corresponding single cluster transition operators. This is done in the next subsection. 
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B. Cluster transition operators 



The partition transition operator Tc can be calculated from the transition operators for the two clusters C\ and 
C2 that constitute the partition C . Two approaches for carrying out calculations of this nature have been presented 
in the literature. First is the approach of Sloan jl5| for a 4-particle system, see his Sec. III. This is essentially 
similar to the method of the last subsection, dividing up the partition transition operator into two parts as in Eq. 
( |l5| ) and obtaining coupled equations for the two parts, as in Eq. (|l8|). The second approach is to separate the 
partition resolvent into cluster resolvents by the use of a convolution. This method is useful since the two clusters 
are dynamically independent, so all operators for one cluster commute with all operators of the other cluster. This 
property of the resolvent was pointed out by Bianchi and Favella p2| and emphasized for use in parts of the 4-particle 
problem by Haberzettl and Sandhas f24\ . Quantities that are used by both methods are defined first, then the methods 
are discussed in turn. 

On the basis that the two clusters in a partition are dynamically independent, the partition potential Vp is a sum 
of cluster components 

V c = V Cl +V C2 . (20) 

Since the kinetic part, H , of the hamiltonian also separates into cluster components, it follows that the partition 
hamiltonian also separates into commuting cluster hamiltonians. That is, these hamiltonians are related according to 
(.7 = 1,2) 

Ho = K Cl + K C2 , H C) = K C . + Y Cj , H c = H Cl +H C2 . (21) 
In a similar manner, cluster transition operators can also be defined according to 

T Cj (z')=V Cj +V Cj -J— T C) (z') 

- V c 3 +Vc,— (22) 

What complex parameter z' is to appear in each cluster transition operator depends on how it is to be used, and 
differs between the two methods. 

In the first method, modelled on Sloan's |]l5| approach, the partition transition operator is written as a sum, 

T c = T Cl + T° 2 , (23) 

whose parts are defined as 

T°' =V Cj +V c .R T c . (24) 



Then repeating a process analogous to deriving Eq. (18), the components of the partition transition operator can be 
shown to satisfy the coupled equations 



\ T r -- ) I T, ,!' ? ! ) ~~ \ T« ,i :•) )**{tC> )' ( 25 ) 
where the parameters 

zi=*-K a , z 2 = z-K 1 (26) 

have been chosen so that the energy factors are consistent with the properties of the partition transition operator. 
According to this, the partition transition operator can be expressed as the series 

T c = T Cl (zi) + T C2 (z 2 ) + T Cl (zi)R T C2 (22) + T Ca (^)R T Cl {z x ) + ■ ■ ■ . (27) 

The cluster transition operators Tcj ( z j ) can be evaluated for arbitrary Zj entirely in the respective cluster subspace 
involving the states of -j or Af ~ 1 particles in Cj . This can be accomplished in the same manner as that described for 
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the N particle system since the cluster is dynamically independent. But it must be evaluated with the appropriate 
Zj parameter when used in calculating the partition transition operator. 

In the second method the basic starting point is to express the partition resolvent as the convolution 

Rc ~2^J r (z'-H Cl )(z-z'-H C2 ) (28) 

of the resolvents of the two clusters. Here the contour T is to be the straight line from — oo to oo, but lying above the 
real axis and below z. Thus z' must be such that < < It follows that the partition transition operator 

is given by 

Tc = V Cl + V C2 + ^ jjy Ci + V C2 ) {z , u Jf_ - - Hca) (V Cl + V Ca ). (29) 

This can be expressed in many different ways. The following emphasizes how the integrand can be expressed in terms 
of the cluster transition operators, Eq. (P4), but leaves the contour integral unchanged. 

The integrand consists of four terms, according to the different cluster potentials. The diagonal in cluster potential 
terms can immediately be recognized as related to the corresponding cluster transition operator, whereas the remainder 
needs the identity 

v ^7^ = T ^7^r (30) 

and its other various combinations. After some calculation, the integrand can be written in the form 

(V Cl +V C2 )— — — — __(V Cl +V C2 ) 

(z' - H Cl )(z - z' - H C J 

1 r _ , 



[T Cl (z') - V Cl ] + [T C2 (z - z') - V C J- 

z — z — rl(7 2 z — tl( 

+Vci (z'-H Cl )(z-z'-U C2 ) Vc2 + z'-H Cl VclVc2 ~- z> - H 

[Tc^O - V Cl ] 1— - + [T C2 (^ - z') - V C2 ] - 1 

z — z- K C2 z' - K Cl 

+[T Cl (^)-Vc 1 ] - ,/ v Tc 2 {z-z')- 



c 2 



>z-z'-K C2 >z-z'-K C2 
' -T Cl (z')T C2 ( Z -z') J——. (31) 



z'-K Cl ' z-z'-Kc. 

Most terms involve the product of the transition operators for the two clusters, weighted with different combinations 
of the free particle resolvents for the clusters. It is up to the computational method to decide which way these are 
to be evaluated. But the first two terms each involves only one cluster transition operator. The contour integral of 
these two terms can easily be done. For the first term the combination Tci ( z ') ~ Vcf I is analytic in the upper half 
z'-plane and vanishes for \z'\ — > oo, so closing the contour at SJ(z') — > +oo contributes a pole only when z' = z — Kp 2 . 
Analogously for the second term, so that the partition transition operator can be written 

T c = T Cl (zi) + Tc 2 (z 2 ) + ^^ {T Cl T C2 terms}. (32) 

While this is only one way of organizing this result there are many other ways in which the cluster transition operator 
expansion of the partition transition operator could be written and no attempt is made here to catalog all the 
possibilities. In comparing the two methods, Eq. ( |2^ ) is an infinite series in cluster transition operators, while Eq. 
( |32| ) is only quadratic in cluster transition operators, but requires an integration over how the energy is divided up 
between the two clusters. 
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In this way the calculation of the partition transition operator has been reduced to the independent computation of 
the cluster transition operators. The advantage of this decomposition is that the individual cluster transition operators 
deal with isolated sets of particles whose number is less than N. Such a cluster of particles can be decomposed into 
partitions as was the original problem, and the whole process repeated. Specifically, for N = 2™ even, the successive 
problems deal with 2 k particles, k = n — 1, n — 2, 1. 



IV. GENERALIZED FADDEEV COMPONENTS 



As a function of z, the resolvent operator R has singularities at the spectrum of the system hamiltonian, thus it 
is of practical use for finding solutions of the Schrodinger equation |29j] , for determining the dynamical evolution of 
the system ]30|-|32]|, and identifying normalizable resonance states [j33[. Faddeev's approach ||| is followed in order 
to express the TV-particle resolvent operator in terms of selected two-cluster partition transition operators. This 
representation of the total resolvent operator is then used to decompose a scattering wave function originating from 
a particular incoming two-cluster partition. The resulting wave functions are here referred to as generalized Faddeev 
components. The homogeneous system of equations associated with the generalized Faddeev components is shown 
explicitly to solve the Schrodinger equation. This section first describes the generalized Faddeev equations for the 
-/V-particle resolvent operator and subsequently discusses the associated wave functions. 

The total resolvent operator R is related to the total transition operator T via the Lippmann-Schwingcr equation 

R = Rq + RqVR = Rq + RqTRq. (33) 



On taking over the partition expansion (15) of T, this can be written as 

R = R + RC > ( 34 ) 

cec 

where the Faddeev type resolvent R c is given by 

R c = R T C R . (35) 
The generalized Faddeev equations ([t8]) lead to the coupled set of equations 

R c = R T c Ro + 

^CsRoTcRoT Ro 

Bee 

= R c R + 2 5a S R T c R s (36) 
Bee 

with partition resolvent operator defined in ( |l2] ) and related to the partition transition operator by 

Rc = Ro + RjVcRc = Ro + RqTcRo- (37) 



For N = 3, Eq. (|3^) is Faddeev's decomposition of the total resolvent operator whose kernel is not connected until 
after one iteration. For N > 3, the kernel of Eq. (|3^) requires one iteration in the same manner providing the kernels 
of the subsystems are already connected. These relations between the transition operators and the resolvent operators 
are analogous to the original Faddeev equations for the three-particle resolvent operator. The set of resolvent equations 
can also be used for the evaluation of statistical mechanical virial coefficients |33],[34| . 

A scattering wave function for the TV-particle system is determined by the total resolvent operator R according to 

| # c ) = lim ieR{E + is) \ <f> c ). (38) 

e— >0 

This is applied here to the selected set C of compatible two-cluster partitions, with | 4>c) a stationary solution of the 
two-cluster partition hamiltonian Hp of energy E, 

H c \<j> c )=E | C ), (39) 

which may be a distorted wave of the two-cluster subsystem. 

The detailed structure of the scattering wave function \J&c) of Eq. ([58]) is now discussed. This is begun by first 
applying the resolvent expansions (|J) and (|3^) to Eq. ([38]), 
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I 4-c) = lim ieTL {E + ie) \ 4> c ) + V lim ieK B (E + ie) \ 4> c ). (40) 

see 

The first term on the right hand side contributes only if C is the channel with all particles free, namely 

lim ieR (E + ie) | C ) =| cf>c)Sc,C N , (41) 

£—►0 

where Cn is the ./V-cluster partition. The second term can be written as a sum of generalized Faddeev components 

lim ieR B (E + ie) \ <j> c ) =\ ipse)- (42) 

e— >0 

According to Eq. (|36|), the generalized Faddeev component | ipse) can be expanded as 

I ipBc) = lim isHb \ <ftc) - hm ieR | 4>c) 

+ V Sa.b'RoTb hm ieR A \ 4> c ) 
Aec 

= I <Pc)Sb,c- I (f>c)dc N ,c + ^2 <5a,bRoT b I IpAC') 

Aec 

= I <Pc)S B ,c- I </>c)<5c„,c + ^a.bRbVb I Vac)- (43) 



Here the identities 



and 



lim ieR B | </>c) = | <Pc)5 B ,c (44) 



R T B = R B V B (45) 



have been used as well as the definition in Eq. ( ] 42[) . Thus the scattering wave function, Eq. (Q), is expressed in 
terms of the generalized Faddeev components (^3[), 

*c) =| 0c)S c ,c N + I M- (46) 
Bee 

For A=3 this expression for the scattering state is exactly the wave function originally derived by Faddeev. Hence 
Eq. (^) is the generalization of Faddeev's scattering state based on a chosen set C of two-cluster partitions of N 
particles. 

It is easily shown that Eq. ( |46| ) satisfies the Schrodinger equation (||). A special case is the homogeneous analog. 
The proof that this formally satisfies the Schrodinger equation is as follows: 

I *c) = I ^ BC > ( 4? ) 
Bee 

BeC \Aec J 

= X fe Sa,b*«Vb [1 + R B V B ] | Vac) J 
sec \Aec / 

= X RoV B [ X ^4,B I ^AC> + I IPBC) ) 

Bee \Aec J 



AC/ 

sec Aec 
= X r oVb I *c) 



Bee 

RoV | * c ). (48) 



Of course a homogeneous solution of the Schrodinger equation occurs only for bound states energies, so such a solution 
plays no role when solving a scattering problem. 

In general, a cluster can be either a bound, or an unbound but interacting, set of particles, with a free particle 
included as a (1-particle) bound state. If both clusters in a partition are bound then the partition is a true asymptotic 
channel. Whether this is or is not the case, the partition wave function has one of two forms, depending on whether 
one of the clusters is or is not a free particle, namely 

|W = (!^' qc> s (49) 

Here | ipc 1 ) is a bound state of the N — 1-particle subsystem and | qc) is the momentum generalized eigenstate 
for the relative motion of the bound and free states, on the basis that the total center of mass momentum of the 
iV-particle system has been removed from discussion. In the second form, the 2nd cluster is now also a bound state, 
which must be explicitly indicated, while qc is again the relative momentum. If a cluster is unbound, but all particles 
interacting, then the corresponding cluster wavefunction must be replaced by ip^ , corresponding to the scattering 
wavefunction from some initial incoming state. It appears intuitively reasonable that this procedure is correct, though 
its complete mathematical justification may be needed since some limits must already have been taken in order to 
define (in general, a product of) scattering states with definite energy as an input, with the next step involving a 
further limit. 



V. THE 4- AND 8-PARTICLE SYSTEMS 



The treatment so far has emphasized how the ./V-particle scattering problem can be broken down into a number of 
problems involving fewer particles. Once a breakdown has been selected, it is then a case of building up the transition 
operator and scattering wavefunction for the .ZV-particle system. The procedure is illustrated by discussing systems 
of 4 and 8 particles. For definiteness it is assumed that there are no bound states for any number of particles. 

For the specific case of there being four particles, these will be labelled 1, 2, 3 and 4. Then according to Eq. (j^) 
there are 3 two-cluster partitions, which are the three pairs of pair particles 12,34, 13,24 and 14,23. Moreover, Eq. 
(^) states that no scaling of the potential is required. Thus the procedure is to first find the two-particle transition 
operator T\2(z) for arbitrary complex z in the upper half plane. On the basis that all particles are the same species, 
this transition operator is the same for any other pair, except for the labelling. The second step is then to find 
the partition transition operator for the two-cluster partition 12,34. This is accomplished according to Sec. IIIB, 
specifically involving either the first approach, Eq. (p7|), or the second approach, Eq. (p9|), with its integrand expressed 
in terms of the pair particle transition operators, Eq. (pl[). Again this calculation needs to be done only once since a 
relabelling immediately gives the partition transition operator for the other two partitions. The last step for getting 
the four-particle transition operator is to solve the 3x3 matrix equations ( |l8| ) and add the results, Eq. (|l5|). Once the 
transition operator is known, the various generalized Faddeev components can be calculated, as described in Sec. IV, 
and the desired scattering amplitude calculated. In this way the four-particle scattering problem is very similar to the 
three-particle scattering problem, involving only three partitions. The extra complexity is in the added structure of 
the partition transition operators. In contrast, the method of Sloan fl5|| involves a 7 x 7 matrix while the Yakubovsky 
[ p"3"[ approach organizes the wavefunction into 18 components. The approach presented here would seem to be both 
simpler and more efficient. 

The 8-particle problem is discussed with the view of clarifying how the proposed approach works for more compli- 
cated systems. Equation (||) states that there are 35 two-cluster partitions, each cluster having 4 particles. Moreover, 
Eq. (|^) states that v& — 1/15, so the potential needs to be scaled by this fraction. This scaling is required when 
finally solving the 35 x 35 matrix, Eq. (fL8|), for the components of the 8-particle transition operator, so enters into 
the calculation of all (the 8-particle, partition, cluster and subpartition and subcluster) transition operators. Now a 
typical 8-particle partition transition operator is needed, which is a combination of two 4-particle cluster transition 
operators. The latter must be found by solving the problem discussed in the previous paragraph, but now with the 
potential scaled by vg. Since there is no further scaling required in solving for the 4-particle transition operator, it is 
only the v$ scaling that must be applied. Thus the detailed procedure that is to be followed is, in sequence: 1) scale 
the pair potential by v$; 2) calculate the 4-particle transition operator for this scaled potential as described in the last 
paragraph; 3) use either of Eqs. ( pTj ) or (£9|) to obtain a two-cluster partition transition operator for the 8-particle 
system, with the second cluster transition operator obtained from the first 4-particle cluster transition operator by 
relabelling; 4) obtain the 34 other partition transition operators by relabelling and solve the generalized Faddeev 
equations (18). 
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VI. DISCUSSION 



A proposal has been made for the generalization of the Faddeev equations 0] to an arbitrary number, N, of particles. 
Essentially this involves selecting a set of partitionings of the N particles into pairs of clusters and expressing the 
scattering properties of the total system in terms of the scattering properties of the clusters. Each cluster is treated 
in the same way so that after successive treatments an individual cluster contains either only a single particle or a 
pair of particles. Although mathematically, other (combinations of) sets of two-cluster partitions lead to a connected 
kernel as well in the present method, the detailed presentation of this paper has stressed the choice of the set of 
two-cluster partitions in which the pair of clusters in any partition is of nearly equal size. This should be the most 
efficient procedure since then there are fewer types of clusters that need to be treated. 

The proposed method requires, in general, a scaling of the potential and expressing the transition operator, and wave 
function, in terms of partition transition operators which are in turn expressed in terms of the transition operators for 
the two individual clusters composing the partition. An essential simplifying feature is that the transition operator 
for an individual cluster can be calculated independently of the presence of the other particles, so is equivalent to 
solving the scattering of a system of fewer particles, which can in general be treated in the same manner as the original 
system. 

For most iV-particle integral equation theories, the starting point is the pair particle interaction or the two-particle 
transition operator, which is then combined in all possible ways with the transition operators of other pairs of particles. 
In order to assure that all particles are connected, that is, no disconnected diagrams appear, it is necessary to iterate 
through a whole sequence of processes. This complicates the description of the A^-particle problem as it appears, for 
example, in the Yakubovsky Jl3[ equations. In contrast, the present method starts by dividing the iV-particle system 
into a selected set of two-cluster partitions. The treatment assures that all clusters are connected to each other. 
There remains the possibility that there is a disconnectionness within a cluster. This is eliminated by repeating the 
procedure for each cluster as if it were a separate scattering system. In this way the formulation has a very simple 
structure which is exactly similar at each stage to the Faddeev equations, in particular reducing to them if N — 3. 

The total resolvent for the TV-particle system is decomposed in a manner analogous to Faddeev's method. This 
decomposition determines the decomposition of the scattering wave functions, viz. Eq. (|4^), into what are here 
called generalized Faddeev components. The sum of the corresponding homogeneous set of generalized Faddeev 
components is explicitly shown to be completely equivalent to the Schrodinger equation. Such a property was explicitly 
demonstrated in Faddeev's theory and Yakubovsky's theory |b|] but has not been explicitly shown in some other 
theories. This property is of particular importance for finding the (discrete) eigenvalues for the A^-particle system and 
for avoiding spurious solutions (for bosons and fermions, an appropriate symmetrization of the Faddeev components 



may be required). In particular, Federbush |3 
systematic study by Glockle and coworkers |27 



first found a spurious solution for the Weinberg equation [g| and a 
|2q ] indicated that most few-body equations admit the existence of 
discrete spurious solutions. However these spurious solutions were not a problem in finding scattering solutions. 

The cluster and partition wavefunctions from all the stages are needed before the total wavefunction can be obtained. 
In general these are nonphysical in that they represent intermediary results for the final calculation. The same is true 
even more so for the Faddeev components. In contrast, if both clusters in a partition represent bound states, then this 
is a valid asymptotic condition and the channel wavefunction does represent a physical state for the corresponding 
scattering system. 
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